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In this Letter, we derive a sufficient condition of synchronizing limit sets (attractors 
and repellers) by using the linear feedback control technique proposed here. There 
examples are included. The numerical simulations and computer graphics show that 
our method work well. 

Historically, the study of synchronization phenomena of dynamical systems has been an active 
topic in physics. In the 17th century, Huggens found two synchronization clocks, other early 
discovered examples, such as, wobbly bridges, the oscillating uniformly Josephson junctions, the 
synchronized lightning fireflies, synchronization of adjacent organ pipes, emerging coherence in 
chemical oscillators, etc [1]. 

In recent decades, the research on synchronization moved to chaotic systems. As we know, 
an essential characteristic of chaotic system is that its evolution sensitively depends on initial 
conditions, intuitionally, this intrinsically defy synchronization. But in 1990 Pecora and Carroll [2] 
pointed out that when we coupled two identical chaotic systems, the synchronization between them 
is possible. In the sequel, lots of methods and techniques of synchronizing two chaotic systems 
were proposed and studied. Amongst these methods, one important method is the linear feedback 
method due to the fact that the drive and response systems become weakly coupled in the process of 
synchronization and this can easily be implemented in circuits. As far as we know, there have been 
no theoretical results available for some interesting chaotic systems to ensure that these systems 
can be synchronized by using the usual linear feedback technique, for example, Rossler system. 
Although it just has one nonlinear term, the Rossler system doesn't have symmetry property 
(Lorenz system has two nonlinear terms, but it has symmetry properties), so it is usually difficult 
to construct a Lyapunov function for proving the global asymptotical stability of the error system. 
In this Letter, we take the Rossler system as an example. The associate synchronization is easily 
realized by utilizing our method. 

On the other hand, when the conditional Lyapunov exponents of two coupled systems are all 
negative, then it is usually thought that these two systems can be synchronized [3]. However, 
it has been recently reported that the negativity of conditional Lyapunov exponents is neither a 
sufficient condition nor a necessary condition for chaos synchronization because of some unstable 
invariant sets in the stable synchronization manifold, see Refs. [4] and [5], and the references cited 
therein. This motivates us to find a sufficient even universal condition of synchronization for more 
chaotic systems. In this Letter, we propose a sufficient condition of synchronization of limit sets. It 
is known that present studies of synchronization are for chaotic attractors (i.e., stable limit sets). 
Unlike current research, in our article, the limit sets are not necessary to be stable limit sets 
(attractors), they even are unstable limit sets (repellers), for example, unstable limit cycle. As 
far as we know, this is the first time to consider synchronization of unstable limit sets. 
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Consider the following system 

where X = (xi, x 2 , ■ ■ ■ , x n ) T 6 R n , / : R™ — » K n is differentiable. 
We build the drive and response systems as follows, respectively, 

^jr = /fai, • • • , Sfe, a;fe+i, ■ • ■ , (2) 

and 

dY 

-^ = f&i, ■■■ ,Xk, J/fe+i, ■ • • , 2/n) + "(y - X), (3) 

in which F = (y 1 , y 2 , ■ ■ ■ , y„) T , F = {x 1: ■ ■ ■ , x k , y k+ i, ■■■ , y n ) T £ R™, < k < n, u is the 
control parameter. Here k = corresponds to one way coupled synchronization approach [6,7]. 

Letting E = Y- X,E = Y- X, and subtracting (2) from (3) yield [8,9] 

= Jo 1 D/(X + s^)^ ds + uE 
Now, we choose a Lyapunov function asF=||i?|| 2 /2, then, 



=< — ,E> = <^Df(X + sE)Eds,E>+u\\E\\ 



dV _ dE 
~di ~ < ~dt 

< \\Df\\-\\E\\-\\E\\+u\\Ef 

< (||A/ll+")l|i?ll 2 , 

where \\Df\\ denotes sup ||D/(X)||, vector norm is 2-norm, matrix norm is the spectral norm. 

So, if ||£>/(X)|| is bounded by a constant M, i.e., ||D/|| < M, then we can choose u < —M 
such that synchronization between systems (2) and (3) can be reached. In more details, if system 

(2) has a (stable or unstable) limit fi, an long as we choose u < —M = — sup then this 

xen 

limit set can be synchronized. Till now, almost all publications are for synchronization of stable 
limit sets, for example, chaos synchronization. But synchronization of unstable limit sets has not 
been studied yet. This paper is the first one to consider such a topic. On the other hand, for 
any continuously differentiable system, any (stable and unstable) limit sets of this system can be 
synchronized with the aid of a simple linear feedback controller derived here. In this sense, our 
derived synchronization method is universal. 

It is evident that the condition of synchronization in this article is sufficient but not necessary. 

Generally speaking, to estimate sup ||D/(X)|| is not easy due to two facts: 1) the bound of the limit 

xeii 

set f2 is often difficult to estimate, 2) the corresponding eigenvalues of (Df(X)) T Df(X) are difficult 
to determine. So we can find a suitable u by numerical exploration such that synchronization can 
be realized. This numerical exploration is described as follows, 1°) choose a suitably large — u such 
that synchronization can be realized, 2°) check u/2 whether or not synchronization can be realized, 
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if does, repeat 2°); otherwise, choose 3m/4 and repeat 2°). This process can not be ended until the 
control parameter | u | is suitable small. Generally speaking, choosing a small | u | is economical. 

We firstly consider synchronization of a unstable limit set. 

It is known that Lorenz system [10] 

dx dy dz 

— = a(y — x), — = rx — y — xz, — — xy ~ bz, (4) 
dt yy h dt dt y ' 

has a unstable limit cycle C (a repeller) if o = 10, b = 8/3, r = 24.5. See Fig. 1. The drive-response 

system is constructed as follows, respectively, 

dx m . , dy m dz rn 

^ — GyVm X m ) : — — fX m y m X m Zrm ^ — X m y m UmZmi 

and 

^ = cr(y s -x s )+u(x s -x m ), ^ = rx s -y s -x s z s + u(y s - y m ), ^ = x s y s -b s z s + u(x s - z m ). 

The simulations are displayed in Fig. 1, where u = — 6 is chosen by numerical exploration. 

For a (chaotic) system, the drive signal can not be randomly chosen, otherwise, the expected 
synchronization can not be implemented. For example, in Lorenz system, if we define variable z 
as the drive signal, the rest as response signals, synchronization between the drive and response 
systems can not be realized by using the usual (Pecora-Carroll) method, for details, see p. 7 of 
Ref. [6]. Here we can realize synchronization by adding a simple controller. The drive-response 
system is constructed below. 

dx m . . dy m dz m 

^VUm ^7n)j ^ Um ^m^rrn ^ ^niVm "m^mi 

and 

^ = v{y s - x s ) + u(x s - x m ), ^ = rx s - y s - x s z m + u(y s - y m ). 

Numerical simulations are shown in Fig. 2. 

Lastly, we study the Rossler system [11], 

dx dy dz 

— = — y — z, —— = x -\- ay. —r = b + z(x — c). 

dt y ' dt y ' dt y ' 

This system is dissipative and has a chaotic attractor when a = b = 0.2, c = 5.7. Rossler 

attractor is somewhat difficult to synchronize by usual methods. Using the present method, this 

chaotic attractor can be easily synchronized. Similar to the above discussion, the drive-response 

configuration is built as follow, 

dx m ^ dy nl dz m . , 

dt dt dt 

and 

^ = ~Vs -z s + u{x s -x m ), ^ =x s +ay s +u(y s -y m ), ^ = b + z s (x s - c) + u(x s - z m ). 

During the process of numerical calculations, we find that chaos synchronization is reached if we 
simply choose u = —6. These simulations are presented in Fig. 3. 
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Figure 1: Synchronization of unstable limit cycle of Lorenz system. Here (a, b, r) — (10, 8/3, 24.5), 
u = —5, the time step-length is 0.001. (a) The unstable Limit cycle C. (b) The evolution of 
synchronization error. Here ei = x s — x m , ei — y s ~ y m , = z s — z m . 
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(b) 

Figure 2: Chaos synchronization of Lorcnz System. Here (a, 6, r) = (10, 8/3, 28), u = —6, the 
time step-length is 0.008. (a) The chaotic attractor. (b) The evolution of synchronization error. 
Here, ei = x s - x m , e 2 = y s - y m - 
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(b) 

Figure 3: Chaos synchronization of Rossler System. Here a = b = 0.2, c = 5.7, u = —6, the time 
step-length is 0.04. (a) The chaotic attractor. (b) The evolution of synchronization error. Here, 

6l — X s Xmi &2 — Us Vm: 63 = Z s 
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